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Abstract 

The influence of magnitoelastic (ME) interaction on the stabilization of long- 
range magnetic order (LMO) in the two-dimensional easy-plane ferromagnetic 
is investigated in this work. The account of ME exchange results in the 
root dispersion law of magnons and appearance of ME gap in the spectra 
of elementary excitations. Such a behavior of the spectra testifies to the 
stabilization of LMO and finite Curie's temperature. 
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I. INTRODUCTION 



It is well known pL H, that the existence of the long-range magnetic order (LMO) in 
2D ferromagnetic is impossible. Formally it means, that in a Heisenberg ferromagnetic the 
integral 

(AM) oc r N(u)kdk, 
Jo 

where N(u) = [exp(u/T) — l] -1 , u> oc k 2 is the law of dispersion of spin waves, which 
determines fluctuations of the magnetic moment, diverges on the lower limit. 

Berezinskii, Kosterlitz and Thouless || |J have shown, that in 2D easy-plane magnet- 
ics at finite temperatures there can exist only a quasilong-range order, that structure is 
determined rather by magnetization vortices, then by spin waves. 

However in || it was shown, that the account of magnetodipole interaction in 2D ferro- 
magnetics results in the root dispersion law of magnons u oc Vk. It causes the convergence 
of the integral (AM) and testifies to stabilization of LMO at T < T c , where T c is the 
temperature of the phase transition. Proceeding from this, it is possible to assume, that 
other relativistic types of interactions can stabilize LMO. For example, in paper |J, it 
is shown, that in an easy-plane 2D antiferromagnetic LMO is stabilized by magnitoelastic 
(ME) interaction. However in this work there is used the Holstein-Primakoff representation 
for spin operators, i.e. quasiclassical representation, which is obviously inadequate to the 
microscopicity of the system in study. 

Therefore it is of interest to investigate a question of stabilization of LMO in 2D easy- 
plane ferromagnetic with the exact account of one-ion anisotropy (OA) and ME interaction. 
Such an account can be made with the use of the Hubbard's operators procedure [[?]. ||. 
This technique allows us to obtain the dispersion equation of hybridized magnitoelastic 
elementary excitations (bound ME waves), valid for arbitrary values of a magnetic ion spin 
(S > 1), temperatures (up to temperature of phase transition), and anisotropy constant. 

II. MODEL SYSTEM 

As a model system, we shall consider an easy-plane 2D ferromagnetic (XOZ-b&sic plane), 
with the Hamiltonian: 

H = -Wl{n- n')Sj nl + § E(^) 2 + A £{(S£) V, + (^) a u„ + 

Z n,n> / n n 

f E 

(S*S z n + S z n S x n )u xz ] + J dv [u 2 xx + u\ z + 2au xx u zz + 2(1 - a)u 2 J, (1) 

where S l n is the spin operator in site n, (3 > is the OA constant, A is the constant of ME 
interaction, Uij is a symmetric part of tensor of deformations, E is the modulus of elasticity, 
a is the Poisson coefficient. To simplify calculations we shall assume, that a magnetic ion 
spin 5 = 1. 

In the Hamiltonian (1) the first two terms describe the magnetic subsystem, the third is 
the two-dimensional ME interaction, and the fourth is the elastic subsystem. 
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As we study the possibility of stabilization of LMO due to the presence of ME exchange let 
us suppose that ME interaction creates the nonzero magnetic moment, which for definiteness, 
we consider parallel to OZ axis. 

After separation of the mean field in the exchange part of (1), we obtain the one- ion 
Hamiltonian: 

H (n) = -I z S z n + ^{Slf + \[u xx {S x n ) 2 + u zz (S z n ) 2 + u xz (S x n S z n + S z n S*)}. (2) 

Solving the one-ion problem with the Hamiltonian (2) we find the energy levels of a magnetic 
ion: 

£i = f + i(4°2 + 2 M (° z ))- x , £ = f + A4°2, £-i = f + i(4°2 + 2^)) + x , 



2\ xx 1 "zz/ A) u 2 xxi — ± 4 1 2 

x 2 = i 2 + (f 4°2 - f ) 2 , /, = /(o) < s' >, 



(3) 



and eigenf unctions of H (n): 

=cos<S|l)+sin<5| -1), *(0) = |0), *(-l) = -sin<5|l) + cos<5| - 1), (4) 

where 

-M^-f)[(x-'.>M^-i)T 1/2 . 

and |i) is the eigenvector of the operator S z , uf) are spontaneous deformations determined 
from the condition of the free energy density minimum: F = F e — T In Z, where F e is the 
density of an elastic energy, Z = Y,m ex v(~E M /T) is the partition function. In our case the 
spontaneous deformations have the form: 

„(°) - . 1 ~ 2a „(°) - . 2 ~ a „(o) _ n 
u xx ~ £j 2 ' ~ E 2 ' ~~ 

Note, that the calculations are carried out in a low temperature limit (the temperature 
is much less than the temperature of the phase transition). In this case it is possible to 
account only the lowest energy level, which as it follows from (3), is E\. 

In the basis of eigenfunctions of the operator Hq(ji) we shall build the Hubbard's oper- 
ators X^ M = \iff n (M'))($ n (M)\, which describe transitions of a magnetic ion from a state 
M' to a state M. These operators are connected with the spin operators by the relations: 

S+ = V2cos 5 ■ (X*° + Xr 1 ) + V2sm5- (X" 1 - X,; 10 ), 

s- = (s+r, (5) 

S z n = cos 25 ■ {Hi - H- 1 ) - sin25 • (H^ - H~ u ), 

where = ± iS%, H^f = X^ M are diagonal Habbard's operators. 

In terms of Habbard's operators the one-ion Hamiltonian can be represented in the form 

H (n) = (EH n ) + AtUdxPtf 1 + K°) + d 2 (X^ 1 + X' 10 )], 
where we denote: 
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d\ = —i=(cos5 — sin 5), d 2 = j={cos5 + sin 5). 

V 2 v 2 

Further, we represent components of tensor of deformations in the form Uij = ufj 1 +u\^\ 

where t£ is the dynamic part of tensor of deformations describing oscillations of a lattice. 

Applying to u[j the method of harmonic quantization ||, from the Hamiltonian (2) we shall 
obtain the Hamiltonian, describing processes of transformations of magnons to phonons and 
vice versa: 



Pi 



Htr — J2ti[J2m PmH^ 1 + J2a ^aX"] ) 

h^ q Ab q ^ + blAT^){q,X). 



(6) 



Here b q! x(bl_ q ^) are rising and lowering operators of phonons with polarization A, T^( a \q, A) 
are the amplitudes of transformations. 

III. SPECTRA OF ELEMENTARY EXCITATIONS 



As it is known, the spectra of elementary excitations of a system are determined by the 
poles of a Green's function, which we shall define as follows: 



G aa (n, r; ri, t" 



■(TX«(t)X«/(t')>, 



where X%(t) = e Hr X^e~ HT , H = H int +H t r, f is a time-ordering operator, and the averaging 
will be carried out with the Hamiltonian H. 

As we work in a mean field approximation, for futher calculations we need only "transver- 
sal" part of the Hamiltonian Hi nt , which in terms of the Habbard's operators has the form: 



n i n' 



n,n';a,/3 



Vector-columns Af and B? have the form: 



27|| (-«) 

r„(M) 

( 7||(«) 
U||(M) 



A a 2 



7l(«) ' 






An 



B< 



'7l(-«)' 




The denominator of the Green's function, which satisfies to the Larkin's equation 
the dispersion equation: 



det 



8a + 



2 -BfY> aa A" + } 



1 - D\(k, u n )TP(k, \)ZW'T-P'(-k, A) 



gives 



0. (7) 



In (7) we denote D\(k,u n ) = 2ui\{k) /[uj 2 — oj\{k)} is the Green's function of noninteracting 
phonons, u)\(k) = C\k is the low of dispersion of A - polarized phonons and c\ is their velosity, 
S aQ (k,u) is the nonreducible (by Larkin) part, T a (k,u) are corresponding amplitudes of 
transformations. In the mean field approximation 
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Z aa '(k,u) = 5 aa ,b(a)G%(L;), 

where G%(u) = [uo + (aE)]' 1 is the zero Green's function, b(a) = (aH) . 

The solutions of the equation (7) determine spectra of quasiphonons and quasimagnons, 
however, as we study the possiblity of stabilization of LMO in a 2D easy-plane ferromagnetic, 
i.e. the question of convergence of the integral determining fluctuations of the magnetic 
moment, first of all we have to obtain spectra of quasimagnons. 

It is necessary to note, that the dispersion equation (7) is valid for arbitrary temperatures 
and values of constant of OA. 

The solution of the dispersion equation in a low temperature limit allows us to determine 
spectra of quasimagnons: 

ep = 2 X , el = [E w - I(k) ■ (1 - sin25)] • [E 10 - I(k) ■ (1 + sin 25)], (8) 

where E w = E\ — E , ep is a high-frequency magnon branch, e a is a low-frequency magnon 
branch. 

We study the solution (8) in the two limiting cases: the case of small and large OA value, 
a). In the case of small OA value, i.e. at j3 <C J , a low-frequency magnon branch has 
the spectrum: 

e 2 a (k) = (b + ak 2 )-(b + ^ + ak 2 ), (9) 

where bo = 3A 2 /4_E is the parameter of ME exchange, a = IqRq, Ro is the radius of 
interaction. In the case of the absence ME exchange we have a spectrum of an easy-plane 
ferromagnetic: 



e a (k) = ak^ + ak 2 . (10) 
b). If OA value is large, i.e. (/3/4) 3> Iq, the spectrum of quasimagnons has the form: 

= (f + «o) • + ao - 2/ + 2ak 2 ). (11) 

Here a® = A 2 (l + a)/2E is the parameter of ME exchange. The spectra of quasiphonons 
are also determined by the equation (7), however, as they do not influence the size of 
the fluctuations of the magnetic moment, they are irrelevant. For example a spectrum of 
longitudinal quasiphonons is as follows: 

uj\k)=u 2 t {k).{l-^). 

The similar result is obtained for t-polarized quasiphonons. It is obvious, that ME interaction 
merely renormalizes the velocity of the sound. 
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IV. THE CRITICAL TEMPERATURE 



Consider fluctuations of the magnetic moment, for example, ((S z ) 2 ). The most simple 
way to calculate it is to represent the operator S z in terms of bose operators through 



bosonisation of the Habbard's operators [10] 



Following [|1(|, we assign to the Habbard's operators X° pseudohubbard's operators X%, 
which are related with rising and lowering bose-operators of quasiparticles by the following 
relations: 

X™ = (1 - a{a n - b{b n )a n , X* = a{, X 1 ' 1 = (1 - a{a n - b{b n )b n , X~ n = b{ , 

yO— 1 f y. V — 10 ht n fjO „t n If~^ h ff^ 1 n h 

a n u n? > tt ni n u n u «' n n u n ni n n 1 a n Uj n u n u n- 

Here a is the bose operator corresponding to the transitions from a state 1 to a state and 
vice versa, b corresponds to the transitions from a state 1 to a state -1 and vice versa. 

Rewriting the Hamiltonian H in terms of bose operators, and using the fact that the 
one-ion Hamiltonian in the basis of its eigenfunctions becomes diagonalize: 

H (n) = (EH n ), 

we write down only two-partical Hamiltonian, that has the form: 



# (2) = J2( E o -Ei- I{k)) ■ a\a k + - E x ) ■ b{b k - \ sin 25 ■ £ I(k)(a k a_ k 



o -- ^ -■ ■ a W-k) 



Diagonalizing the obtained Hamiltonian by the standard u — v transformation: 

a\ = u k a k + v k a ] _ k , a k = u k a\ + v k a- k , 

we receive: 

H {2) = ^e a a{a k +J2 £ /3^k, 

k k 

where expressions for bosone energies coincide with (8). 

Using the connection of spin operators with the Habbard's operators (5), for fluctuations 
of the magnetic moment we shall receive: 

„„.,, ..,+ , . 1 — sin25 , + . l + sin25 . , 

((^) 2 > = - sin 28(b%) + (4 a n ) + . (12) 

Summarizing (12) we receive the desired average and 

1 \~^/iA'L \ ^ [°° kdh 1 / t \ ^ [°° kdk 

where e a (k) and £p{k) are determined by the expressions (8). The relation (12) is valid for 
arbitrary values of OA constant. As in the previous case we investigate two limiting cases: 
small and larg OA values. 
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a). In the case of small OA the contribution of a high-frequency branch of excitations to 
fluctuations of the magnetic moment can be neglected. It is obvious, that in this case the 
integral 

kdk 



converges on the lower limit (e a (k) is determined by the formula (9)), hence, in this case 
there exist a LMO in 2D ferromagnetics. If we disregard with the ME exchang, the integral 
diverges on the lower limit, that testifies to the absence of LMO. Besides from the condition 
(1/N) J2 n ((Sn) 2 ) = 1 it is possible to determine the temperature of the phase transition. In 
our case we have: 

kdk 



(2vr) 2 Jo e £ «/ T - - 1 

where T c is the temperature of the phase transition. Substituting the expression (9) for 
e a (k), we obtain the expression for T c : 



In 



vV(&o + /?/2) 



(13) 



As it is evident from (13), the temperature of the phase transition is determined both by the 
ME interaction, and OA. But still ME interaction is a decisive one, and at b = 0, T c — > 0. 
b). In the case of large OA value fluctuations of the magnetic moment is equal to 

1 \^/(qx\2\ __ I f 00 kdk 
JVV (27r) 2 7o e WT_i' 

where e a (k) is determined by the formula (11). Bisides even in the absence of the ME 
exchange the integral does not diverge on the lower limit (in a spectrum of quasimagnons 
there is a finite gap at a = 0). Evidently, such a behavior of fluctuations is connected with 
the fact that in the case of large OA there may realize the phase with a tensorial order 
parameter ]TIJ (the so-called QU-phase). 

As before, we shall determine T c . In our case we have (/3/4 ^> I ): 

Te^^, Zna.p, A«^-2/ + ao). (14) 



V. CONCLUSION 

The carried out investigations show, that the account of the ME interaction results in 
the stabilization of the LMO in a 2D ferromagnetic. From the formula (13) it follows, that 
in a low-anisotropic FM the value of T c is determined by the ME gap of a quasimagnon 
spectrum. 
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In the case of large OA the influence of ME interaction on the establishing of LMO is 
not so essential, since in such systems there can exist a magnetically ordered state with a 
tensorial order parameter (the so-called quadrupole phase) ||12|| . Temperature of transition, 
as follows from (14) is determined mainly by the constant of OA, and the ME interaction 
only renormalizes it. 

One of the authors (S.D.) thanks International Soros Program of support of education 
in the field of exact sciences, Grant Nr SU072163. 
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